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We study the effects of Majorana modes and interactions between electrons on transport in a 
one-dimensional system with a junction of three p- wave superconductors (SCs) which are connected 
to normal metal leads. For sufficiently long SCs, there are zero energy Majorana modes at the 
junctions between the SCs and the leads, and, depending on the signs of the p-wave pairings in the 
three SCs, there can also be one or three Majorana modes at the junction of the three SCs. We show 
that the various sub-gap conductances have peaks occurring at the energies of all these modes; we 
therefore get a rich pattern of conductance peaks. Next, we use a renormalization group approach 
to study the scattering matrix of the system at energies far from the SC gap. The fixed points of 
the renormalization group flows and their stabilities are studied; we find that the scattering matrix 
at the stable fixed point is highly symmetric even when the microscopic scattering matrix and the 
interaction strengths are not symmetric. We discuss the implications of this for the conductances. 

Finally we propose an experimental realization of this system which can produce different signs of 
the p-wave pairings in the different SCs. 


1 Introduction 

The subject of Majorana fermions has been exten¬ 
sively studied in recent years, both from the fundamental 
physics point of viewi as well as for possible applications 
in topological quantum computation^. A well-studied ex¬ 
ample of a system which hosts Majorana modes is the 
Kitaev chain, which is a spin-polarized p -wave supercon¬ 
ductor (SC) in one dimension. This system is known to 
have topological phases in which there is a zero energy 
Majorana mode at each end of a long chain^. This system 
and others similar to it have been theoretically studied in 
a number of papersi^~—. Several experiments inspired 
by these models have looked for Majorana modes, either 
through a zero bias conductance peak^r— or through the 
fractional Josephson effect^. The zero bias peak occurs 
because the Majorana mode (which lies at zero energy 
for a long enough system) allows an electron to tunnel 
from a normal metal lead into the superconductor where 
it turns into a Cooper pair; this gives rise to a perfect 
Andreev reflection back into the lead. 

The Kitaev chain has been generalized in a number 
of ways; some of these generalizations give rise to more 
than one Majorana mode at each end of the SC re- 
gior&^ i 22 ! 37 ! 41 . The case of additional Majorana modes 
appearing inside the bulk of the SC region (rather than at 
the ends) is less well studied. Strictly speaking, the term 
Majorana modes is used to refer to states with exactly 
zero energy. In this paper, we will discuss various sub¬ 
gap modes which lie within the SC gap but may or may 
not lie at zero energy. All these modes originate from the 
zero energy Majorana modes which occur if the SC wires 
are very long, but they hybridize and thereby move away 
from zero energy when the wire lengths are finite. It is 
interesting to study the effect of all such sub-gap modes 
on the electronic transport across the system. 

The effect of interactions between the electrons is also 
of interest. In one dimension, interactions are known 


to have a dramatic effect, turning the system into a 
Tomonaga-Luttinger liquid. The fate of the Majorana 
modes in the presence of interactions has been studied 
extensivel^i^^^^d^^ We will examine a dif¬ 
ferent aspect of our system, namely, the fixed points of 
the conductances which arise from the renormalization 
group flows induced by interactions. 

In this paper, we will study the conductances of a 
one-dimensional system consisting of a junction of three 
finite-sized p-wave superconductors, each of which is in 
turn connected to a semi-infinite normal metal (NM) 
lead. We will study transport in two different regimes 
of the energy, namely, inside the SC gap and far outside 
the gap. 

In the first part of this paper, we will study the Ma¬ 
jorana and related sub-gap modes and their effect on 
the sub-gap conductances. A motivation for studying a 
junction of three p-wave superconductors is that a topo¬ 
logical quantum computation using Majorana modes in 
quasi-one-dimensional systems would involve braiding of 
two Majorana modes without the two modes ever com¬ 
ing close to each other. In order to achieve this we must 
consider junctions of three or more wires, so that a Ma¬ 
jorana mode can be moved from the first wire to the sec¬ 
ond without disturbing a Majorana mode lying on the 
third wire. (A single superconducting wire does not al¬ 
low braiding of Majorana modes). We will show that 
there are two interesting possibilities for a junction of 
three superconducting wires which have rather different 
properties: the p-wave pairing amplitude A can have the 
same sign in all three SCs, or A can have the same sign 
in two of the SCs and the opposite sign in the third SC. 
We will see that the number of zero energy Majorana 
modes is different in the two cases when the SC wires are 
very long (namely, much longer than the decay length 
which will be discussed below). In the first case, there 
are six Majorana modes (three at the NM-SC junctions, 
and three more at the junction of the three SCs), while in 
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the second case, there are four Majorana modes (three at 
the NM-SC junctions and only one at the junction of the 
SCs). The fact that there can be three Majorana modes 
at a junction of three p -wave superconducting wires has 
not been pointed out earlier as far as we know. 

We would like to note that many kinds of junctions 
have been studied before, such as junctions of three 
wires or three spin- 1/2 chains (which can be mapped to 
fermionic chains) 11 ! 66 ’ 67 , a superconducting junction of 
non-superconducting wire o 68 ’ 69 , and a junction of mul¬ 
tiple non-superconducting wires with a superconducting 
wire 70 . In this paper, we will not consider the effect of 
interactions at the junction (i.e., charging energy) which 
leads to the Kondo-Majorana physics studied in Refs. [ 68 | 
andjU. However we will discuss the effect of interactions 
in the wires away from the junction as mentioned below. 

Electronic transport across a NM-SC-NM system has 
been studied for many yearsZIr— . The presence of a 
SC means that there will be both normal reflection 
and transmission and Andreev reflection and transmis¬ 
sion^^. Hence there are two kinds of differential con¬ 
ductances which can be measured in this NSN system: 
a conductance from one NM lead to the other NM lead 
which we will call Gjv, and a Cooper pair conductance 
from a NM lead to the SC which we will call Gc■ In our 
system with three SCs and three NM leads, we will con¬ 
sider an electron incident from one of the leads (called 
NM1); we can then have a conductance Gc from that 
lead to the SC, and conductances Gjy 2 and Gpj 3 from 
that lead to the other two leads called NM2 and NM3 
respectively. We will use a continuum model for this sys¬ 
tem, rather than the Kitaev model which is defined on a 
lattice. We will first present the boundary conditions at 
the three NM-SC junctions which follow from the conser¬ 
vation of both the probability and charge currents. We 
will then discuss the boundary conditions at the junc¬ 
tion of the three SCs; this will turn out to involve a 3 x 3 
Hermitian matrix M which determines how a current in¬ 
cident on the junction from one of the SCs either gets re¬ 
flected back to that SC or gets transmitted to the other 
two SCs. Using all these boundary conditions, we will 
numerically calculate Gc, Gn 2 and Gn 3 as functions of 
the energy E of the electron incident from the lead NM1 
and the lengths Lj of the three SCs. We will see that 
Gc has a rich structure of peaks when E lies in the su¬ 
perconducting gap. The conductance calculation will be 
followed by the discussion of a box made of only the SCs 
with hard wall boundary conditions (namely, with no NM 
leads). We will numerically calculate the energies of the 
sub-gap modes in this system and show that this explains 
the locations of the peaks in the conductance Gc of the 
system with NM leads. 

In the second part of the paper, we will study the effect 
of interactions between the electrons on the conductances 
of the system at energies which are far from the SC gap, 
namely, when \E\ |A|. (The Majorana and other sub¬ 

gap modes play no role in this part; this is therefore 
complementary to the first part of the paper where we 


look at the sub-gap conductances). It is common to use 
the technique of bosonization to study one-dimensional 
systems with interacting electrons^!. However, for rea¬ 
sons that will be discussed below, it turns out to be diffi¬ 
cult to use bosonization when there is a junction of three 
or more wires and superconductivity is present. (Three- 
wire junctions without superconductivity have been stud¬ 
ied in Refs. [7814821 using bosonization and in Ref. [H 
using functional renormalization group methods. Junc¬ 
tions between a superconducting wire and multiple non¬ 
superconducting wires with interactions have been stud¬ 
ied using bosonizationSi). We will therefore use a dif¬ 
ferent method which is valid when the strength of the 
interactions is weah2£~—. Using this method we will find 
that the scattering matrix S which characterizes the junc¬ 
tion of three wires effectively becomes a function of the 
length scale, and we will find a renormalization group 
(RG) equation for S. Using the RG equation, we will 
study how S varies with the length scale; we will find the 
fixed points of the RG equation and study their stabili¬ 
ties. We will then discuss the implications of this for the 
conductances of the system when the wire lengths are 
large or the temperature is low. 

The plan of the paper is as follows. In Sect. 2, we 
introduce the model for the NSN system and derive the 
boundary conditions at the junctions between the SCs 
and the NM leads and at the junction of three SCs. We 
show how this can be used to derive the various differ¬ 
ential conductances Gnj and Gc at energies lying inside 
the SC gap. In Sect. 3, we numerically calculate G^j and 
Gc for two cases: when the p- wave pairings A j have the 
same sign in all the three wires and when one of them 
has a different sign from the other two. We discuss three 
different regimes of the lengths of the SC wires: these 
lengths can be less than, a little larger than, and much 
larger than the decay length of the sub-gap modes which 
appear near the different junctions. The locations of the 
peaks in the conductances are quite different in the three 
length regimes and also in the two cases of the relative 
signs of the A j. To understand these differences, we have 
numerically found the energies of the sub-gap modes and 
shown that they precisely match the locations of the con¬ 
ductance peaks. In Sect. 4, we provide analytical argu¬ 
ments to show that the number of zero energy Majorana 
modes at the junction of three long p- wave SCs is three 
if the Aj’s have the same sign but is one if one of the 
Aj’s has a different sign. In Sect. 5, we study the effect 
of interactions between the electrons on the conductance 
at energies far from the SC gap (but much smaller than 
the band width of the system). To this end we consider 
a junction of three NM wires which meet in a SC region, 
and we derive the RG equations for the scattering matrix 
of this system for the case where the electrons interact 
weakly with each other. In Sect. 6 , we study the fixed 
points and stabilities of the RG equations and show that 
this can lead to power-law dependences of the conduc¬ 
tances on the wire lengths or temperature. In Sect. 7 we 
show how a system of three p -wave SCs can be experi- 



3 


mentally realized and how the two cases of all the Aj-’s 
having the same sign or one of them having a different 
sign can be fabricated. We end in Sect. 8 with a summary 
of our results. 


everywhere; wc will set Aj = 0 in the NM leads. (We will 
generally set h = 1 in this paper, except in places where 
it is required for clarity). The Heisenberg equations of 
motion idtc = — [H, c] and idtd = — [H, d\ imply that 


2 Model for a system of three SC wires 
making a W-junction 

We begin with a continuum model for a Y -junction 
of three p -wave SC wires in one dimension as shown in 
Fig- HI A NM lead is attached to the end of each SC where 
there is a barrier modeled by a (5-function potential. Each 
NM-SC system has a coordinate system Xj, with j = 
1,2,3; we define the point where the three SCs meet as 
the origin Xj = 0 for all j. As we move away from this 
junction, Xj will be taken to increase. The different SCs 
may have different lengths and lie in the regions 0 < Xj < 
Lj\ the leads lie in the region Lj < Xj < oo. 



FIG. 1: Schematic picture of a T-junction of three wires la¬ 
beled as 1-3. The inner (lighter) regions of the wires are 
p-wave SCs while the outer (darker) regions are NM leads. 
For an electron incident from NM1, six amplitudes are shown 
in the leads: r n i, r a l are normal and Andreev reflections in 
NM1, t n 2 , t a 2 are normal and Andreev transmissions in NM2, 
and t.„ 3 , ta .3 are normal and Andreev transmissions in NM3. 


Let us denote the wave function in wire j as ip = 
(c, d) T , where c(xj,t ), d(xj,t ) are the electron and hole 
components respectively (we assume that the electrons 
are spin-polarized and will therefore ignore the spin la¬ 
bel). The Hamiltonian in each wire can be written as 


i3 “ = -(IC") c “ T7 d *‘ d - 

m = (IC ") d - '-t 3 -> c - <2) 

For a wave function which varies in space as e ±lkXj , the 
energy is given by ±[fc 2 /( 2 m) — p] if Aj = 0 , and by 

±y /, [fc 2 /(2m) — p] 2 + A 2 {k/kp ) 2 if A j 0. The corre¬ 
sponding wave functions will be presented below. We see 
that the energy spectrum in the j-th SC has a gap equal 
to 2 A j at k = ±kp. 

Let us define the particle density p p = c^c + d^d (this 
counts electrons and holes with the same sign) and the 
charge density p c = c^c — did (which counts electrons 
and holes with opposite signs; we are ignoring a factor of 
electron charge here). Using Eqs. J2|) and the equations 
of continuity dtp v + d x J P = 0 and dtp c + d x J c = 0, we 
find the particle and charge currents to b o 71 ' 88 

J v = -f— [— c^d x c + dxC^c + d)d x d — d x d)d\ 

2m 

+ -— {c^d + die), 

Kf 

J c = Ji + [ dx’ Ji(x'), 

Jo 

Ji = J— [~ ( ld x c + d x c*c - dld x d + d x dld\, 

2m 

2A 

J2 = - 7 — (dxdd + dldxc). ( 3 ) 

kp 

The last term, J 2 , can be interpreted as the contribu¬ 
tion of Cooper pairs to the charge current; note that it 
vanishes in the NM where A = 0. 

The boundary conditions at the NM-SC junctions at 
Xj = Lj can be found by demanding that the currents J p 
and J c be conserved at those points. At the junction 
X\ = L 1 , let us consider the wave functions ip nm i = 
(cnmi, d nm i) T and ip scl = (. c scl , d scl ) T at the points 
x\ = L 1 + e and x\ = L\ — e, i.e., in the NM and SC 
regions respectively. The condition J P (L 1 + e) = J P (L 1 — 
e) implies that 


r°° h 2 d x , h 2 d 2 

H= dxj [ c t(-- — p)c — d\ -- — p)d 

J 0 3 1 2m J y 2 m 

~ l -^{Jd xl d+dld Xj c)\, ( 1 ) 

where p is the chemical potential, kp = y/2mp/H is the 
Fermi wave number, and Aj is the p-wave superconduct¬ 
ing pairing amplitude which will be assumed to be real 


[ C nmld X i c nml + d Xl C nml C nrn i 

A d nm ^d Xl d nrn 1 d Xl d nrn ^d nm i\ 

[ C sc id xi c sc i 0 Xl c sc ^c sc 1 

d“ d scl d Xl d sc i d Xl dl cl d sc i] 
( c Li^sci + d\ cl c sc i). 


( 4 ) 
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The simplest way of satisfying this condition is to set 


Cscl 

dsci 

f) I , 

Ox\Cscl I ^scl 

Vf 

a , _ *Ai 

0^X1^501 ^-Scl 

V F 


Cnm 1 , 
dnmli 

&x\ Cnml 5 


(5) 


where vp = k F /m is the Fermi velocity. The first two 
equations above mean that the wave function is contin¬ 
uous while the last two equations imply that the first 
derivative is discontinuous in a particular way. We now 
find that Eqs. © imply that the charge current is con¬ 
served, namely, J C {L\ — e) = J C (L i + e). Next, let us con¬ 
sider what happens if a (5-function potential of strength 
A is present at the junction at x\ = L \; note that the 
dimension of A is energy times length. (This potential 
is physically motivated by the fact that in many exper¬ 
iments, the NM leads are weakly coupled, by a tunnel 
barrier, to the SC. This can be modeled by placing a 5- 
function potential with a large strength at the junction). 
Now there will be an additional discontinuity in the first 
derivative at Xi = L i; this is found by integrating over 
the (5-function which gives 


^xi^nml(-blTe) (9xiVbcl(-bl e) — 2(aA (0). (6) 

Hence Eqs. © must be modified to 


C-scl 

dsci 


C-nml 7 
^nml 7 


&xi C-scl 
&x\d S cl 


zAi 

v F 

iAi 

v F 


dscl + 2?77/A Cscl — 

Csci + 2mA d sc i = d Xl d nrn \. (7) 


For the other two NM-SC junctions at x<i = L 2 and X 3 = 
L 3 , we will get boundary conditions similar to Eqs. ©. 

We have found the boundary conditions at the NM-SC 
junctions. Now we look at the junction Xj = 0 where the 
three SCs meet. We can find the boundary conditions 
at this junction using the conservation of the particle 
current J p \ this implies that 

F Jpi = 2^ F [~ C l c jdxn c scj + d x .cl cj c scj ] 

o j 

^ 2 rn FI [dlcidxjdscj - 9 Xj d\ c -d sc j\ 

3 

+ ( C tcjdscj + d\ C jC SC j) 

= 0 . ( 8 ) 


The simplest way of satisfying this condition is to set 

C^xCsc — * Csc T N • dsc, 

d x d sc = M* ■ d sc + N* • C sc , (9) 


where c sc — {csci, Csc 2 , Csc 3 ) , d S c — (d sc i, d SC 2 , d SC 3 ) , 
and M, N are 3x3 matrices. Substituting the above 
conditions in Eq. ©, we find that we must have 


N 


M f = M, 
nT = _i 2 m ^ 
Kf 


( 10 ) 


/ Ar 0 0 \ 

where A = 0 A 2 0 .In our calculations we will 

V 0 0 A 3 / 

assume that the magnitudes of all the A j ’s are the same 
but the signs of A j may vary with j. Now we will use 
the conservation of the charge current J c at the junction 
Xj = 0. The second term in the expression in Eq. © for 
the charge current, namely, j* 3 dx' J 2 (x') goes to zero 
as Xj —> 0. So the conservation of J c implies that 


F — 2m Ft C L jdxjCscj + d Xj c\ c ^c S cj} 

3 3 

+ ^ F i~dicjd Xj d SC j + d Xj d\ C3 d SC j] 

3 

= 0 . ( 11 ) 

Using Eq. © in the above expression we find that 


N = — ^ A. (12) 

k _f 

Hence Eq. © becomes 

ZTfl 

^x^sc — -M- * Csc 1 ^ ‘ d sc , 

Kp 

ZTTl 

d x d sc = M* • d sc + A • c sc - (13) 
k F 

It is useful to note some symmetries of our system. 

(i) Eqs. © are symmetric under time reversal (which 
changes t —> — t and complex conjugates all numbers) 
if we transform c —> c* and d —d*. This will also 
be a symmetry of Eqs. (1131) if M* = M. Eq. (1101) then 
implies that M must be both real and symmetric. We 
will assume this henceforth. 

(ii) Eqs. 0. © and m remain invariant under the 
transformation 


c(xj) —> i c(xj), d(xj) —> —id(xj), 

Aj - A j, (14) 

for all value of j and Xj. This implies that the conduc¬ 
tances discussed below remain invariant if A.,- —> —Aj on 
all the wires. 

We now use the boundary conditions discussed above 
to find the various reflection and transmission amplitudes 
when an electron is incident from, say, the NM1 lead with 
unit amplitude. In the presence of the SCs, the various 
scattering processes that may occur are as follows^. 

(i) an electron can be reflected back to the NM1 lead 
with amplitude r n \. 
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(ii) a hole can be reflected back to the NM1 lead with 
amplitude r a Charge conservation then implies that a 
Cooper pair must be produced inside the region SCI. 

(iii) an electron can be transmitted to the NM2 lead with 
amplitude t n 2 . 

(iv) a hole can be transmitted to the NM2 lead with 
amplitude t a 2 - (This is usually called crossed Andreev 
reflection^). Then charge conservation implies that a 
Cooper pair must be produced inside SC2. 

(v) an electron can be transmitted to the NM3 lead with 
amplitude t n 3 . 

(vi) a hole can be transmitted to the NM3 lead with 
amplitude t a 3 which implies that a Cooper pair must be 
produced inside SC3. 

If the energy E of the electron (incident from the lead 
NM1) lies in the superconducting gap, i.e., —A < E < A 
{E can be interpreted as the bias between the chemical 
potentials of NM1 and the SCs), Eqs. 0 imply that the 
wave functions in the different regions must be of the 
form 

Vw = e-^ kF+k ^ ( l ) + r„i e^ kF+k ^ ( J ) 

+ r a \ e<(-**+*)*i f ° ^ , 

1pnm2 = t n 2 e^ +k ^ ( J ) + i„ 2 e**"*'^ ( J ) , 

VW = t n3 E {kF+k ^ ( J ) + t a3 e*(-^+fc> 3 

^ scj = tlj ^ 3 ( sgn(Aj) e* ) 

+ hj ( sgniAj) e-* ) 

+ ^ ^ ( sgn(Aj) e-** ) 

+ * ^ k4Xj ( -s 5 n(i,) ) ’ < 15 > 

where = (E — i\J A 2 — E 2 )/A, and sgn denotes the 
sign function. (We have assumed that the magnitudes of 
the A j are the same; hence e 1 ^ is the same for the three 
SCs). The top and bottom entries in the wave functions 
denote the particle and hole components. The wave func¬ 
tions in the NM leads are proportional to e l ( ±k F+ k ) x i . 
where k -C kp, namely, we are working close to the Fermi 
energy. In each SC, we have four modes; two of these de¬ 
cay exponentially while the other two grow as we move 
away from the junction of the three SCs. We denote the 
wave numbers of these modes by fci, — k%, ks and — k^. 
Defining the decay length 

f = - rF (16) 

Ay/1 — (E/A) 2 

we find that the decaying modes have 

hi = kF + i/£, 

~k 2 = -k F +i/£, (17) 


while the growing modes have 

k 3 = k F - i/£, 

—k 3 = —k F — i/£. (18) 

From Eqs. 0 we get the relation between the wave 
functions on the NM and SC sides at each of the 
NM-SC junctions. Similarly, from Eqs. m and m 
we get the relation among the wave functions of the 
three SCs at the junction where they meet. We 
thus have eighteen equations for the eighteen unknowns 
Ail? Dzi? t 7 t a 2 , t n 3 , t a 3 , t\j , t 3 j and ^ 4 ^, where 
j = 1,2,3. After solving these equations we can calcu¬ 
late the reflection and transmission probabilities. The 
conservation law for the probability current implies that 

\r n l\ 2 + |r Q l| 2 + \tn 2 1 2 + \ta 2 \ 2 + |^n.31 2 + |^a.31 2 = 1- 

(19) 

The net probabilities for an electron to be transmitted 
from the NM1 lead to the NM2 and NM3 leads gives 
the differential conductances Gn 2 and Gjv 3 respectively, 
where 

Gn 2 = \tn2 1 2 — \ta2\ 2 , 

Gn 3 = \tn3 1 2 — l^a.31 2 - (20) 

The net probability for the electron to be reflected back 
to the NM1 lead is 

Gb = M 2 - |r Q i| 2 . (21) 

The remainder, denoted by the differential conductance 
Gc, is the probability for the electron to be transmitted 
into the SCs in the form of Cooper pairs. The conserva¬ 
tion of charge current implies that 

Gc = 1 — Gjv 2 — Gn3 — Gb 

= 2 (M 2 + M 2 + M 2 ), (22) 

where we have used Eqs. (HUD to derive the last line 
in Eq. ([22]) . [Actually, the differential conductances into 
the NM leads 2 and 3 and into the SCs are given by 
er I (2 -kK) times Gjv 2 , Gn 3 and Gc respectively, where e 
is the charge of an electron. However, we will ignore the 
factors of e 2 /{2 'kK) in this paper and simply refer to Gnj 
and Gc as the differential conductances.] 

We note that a differential conductance denotes G = 
dl/dV. To measure Gn 2 , Gn 3 and Gc in our system, 
we have to assume that there is a voltage bias V between 
the NM1 lead on the one hand and the SCs and the NM 
leads on the other (the SCs, NM2 and NM3 are taken to 
be at the same potential). Namely, we choose the mid¬ 
gap energy in the SCs as zero, and the Fermi energies 
in the NM1 lead as E = eV and in NM2 and NM3 as 
zero. The differential conductances Gnj = dlNj/dV and 
Gc = dlc/dV are then the derivatives with respect to 
V of the currents measured in the NM leads and in the 
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SCs. 

3 Numerical results 

3.1 All Aj’s with the same sign 

In this section we present numerical results for Gc and 
Gn 2 as functions of the length L\ of SCI and the ratio 
E/A lying in the range [—1,1]. The length scale associ¬ 
ated with the SC gap is r/ = vf/A = fcp/mA. (This is 
different from the length £ introduced in Eq. (1161) which 
depends on the energy E. Note that £ = rj if E = 0). We 
study three cases, namely, L < r), L ~ rj and L > rj. 

The values of the parameters that we have used to 
numerically calculate the conductances are as follows: 
kF = 1, m = 0.5, A = 5, and Ai = A 2 = A 3 = A = 0.1. 
Throughout this section and in the next, we will take the 
matrix at the junction of the three SC wires to be of a 
form which is completely symmetric under any permuta¬ 
tion of the three wires, 


The diagonal terms of M connect d x c (d x d) and c ( d ) 
in the same wire, while the off-diagonal terms connect 
d x c id x d) and c (d) in different wires. If we choose the 
off-diagonal and diagonal elements of M to be different, 
some of the results for the conductances may differ from 
what we will discuss below. 

In the first row of Fig. [2] we have shown top views of 
Gc- In the third row we show surface plots of Gn 2 which 
is the conductance in NM2. In our calculations we have 
chosen L 2 and L 3 not to differ much from each other, so 
that Gn 2 and Gn 3 are similar. Hence we have presented 
G N 2 only. 

To understand the conductance plots better we con¬ 
sider another system which is made of three SC wires 
but is not connected to any NM leads; namely, there is 
an infinite barrier at the ends of the SC wires. Thus this 
system is similar to a particle in a box but the box is 
made of three SC wires. The wave function is zero at the 
ends of the wires because of the infinite barriers present 
there. At the junction of the three wires we use the same 
boundary condition that we derived in Eqs. m- Using 
the boundary conditions we get a set of twelve linear ho¬ 
mogeneous equations for the amplitudes t\j , t%j , t$j and 
tij in the SC wires. These equations have a non-trivial 
solution if the determinant D of the matrix constructed 
from these equations is zero. We plot the determinant D 
as a function of L\ and E/A. The points where D = 0 
give the parameter values where sub-gap modes appear 
in the system. In the second row of Fig.[2j we have shown 
the top views of D; the lightest regions correspond to the 
points closest to D = 0. We expect that the positions of 
the peaks of Gc will match the positions of the zeros of 
D as Li and E/A are varied. 

For the parameter values we have chosen, we find that 


?7 = 20. In the first column of Fig. [2] we take L\ = 3.57T 
to 7.5 tt (so that L\ < rj), L 2 = 6.77r and L 3 = 6.37T. 
In this regime we find that there are six sub-gap modes 
inside the SCs which is evident from the top view of Gc in 
Fig. [2] (a). If we look at a particular value of L\ in that 
figure, we see six different modes at different energies. 
(Some of the modes look fainter than the others). Three 
of the sub-gap modes lie near the NM-SC junctions as we 
know from earlier papers (see, for example, Ref. [55l l . The 
other three modes must lie near the three-wire junction. 
To explain the presence of modes at non-zero energies, 
we recall that the momenta in the SC have an imaginary 
part, giving rise to factors of ztxj// as we can see from 
the expression for ip SC j in Eq. (Tl5l) . On wire j, the sub¬ 
gap wave functions of the form e Xj A decay exponentially 
as we go away from the NM-SC junction, while the wave 
functions of the form e~ Xj A decay as we go away from 
the junction of three SC wires. When the lengths of the 
SCs are small compared to the decay length £ of the 
sub-gap modes, the amount of decay of the modes inside 
the SC will be small. Hence, the sub-gap modes at the 
NM-SC junctions and at the junction of three SC wires 
will hybridize and their energies will split from zero. So 
the mixing of the sub-gap modes for SC wires with small 
lengths is responsible for the appearance of states at non¬ 
zero energies. These energies oscillate with L\ and vanish 
at certain values of L\. 

As there are states inside the SCs and the lengths of 
the SCs are smaller than the decay length £, there are two 
possibilities for an electron coming from the lead NM1. 

(i) It can enter the SCI by coupling to the Majorana 
mode sitting there and then turn into a Cooper pair; 
a hole goes back into NM1 to conserve charge. In this 
process we get a finite r a 1 and hence a finite Gc- 

(ii) It can enter the SCI similarly as described above 
and can then get transmitted to the other two NMs as 
the wave functions of the sub-gap modes decay very little 
inside the SCs. So we get finite transmission probabilities 
and hence a finite Gn 2 and Gn 3 - From the figures in the 
first column of Fig. (2j it is clear that both Gc and Gjv 2 
are appreciable when L\ < 77. 

In the third column of Fig. [2J we choose L\ = 70.57T to 
74 .57T, L 2 = 68.37T and L 3 = 68.77r, so that all the lengths 
are much larger than the decay length £ of the sub-gap 
modes. In this regime the sub-gap modes are mostly lo¬ 
calized at the ends of the SCs and at the junction of the 
three SC wires. So the coupling between the different 
sub-gap modes are small, and their energies remain al¬ 
most at zero. In Figs. [2] (c) and (i), we see that instead 
of six sub-gap modes at different energies, we now have 
all the sub-gap modes near zero energy (namely, they are 
almost Majorana modes). There is a small splitting in 
the energy around E = 0 because of a small non-zero 
coupling between the different modes. If we increase the 
various lengths more, we expect to see no splitting at all 
and all the modes should stay exactly at E = 0. 

In the large length regime we find that Gc almost ap¬ 
proaches its highest value of 2 (in units of e 2 /h) while 
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Gn 2 is very small. The reason for this is that the sub¬ 
gap modes are now almost decoupled from each other; 
due to the exponential decay of their wave functions, the 
probability of an electron to travel through the SCs and 


get transmitted to the two NMs on the other side is very 
small. So the electron mostly transmits into the SCs and 
turns into a Cooper pair inside the SCI, so that |r a | 2 ~ 1. 



FIG. 2: Plots of conductances and energies of sub-gap modes when A, has the same sign in the three SC wires: the parameters 
chosen are kr = 1. m = 0.5 and Ai = A 2 = A 3 = A = 0.1, so that the length scale 17 = 20. In figures (a-c) and (g-i), the 
strength of the NM-SC barriers is A = 5, while in figures (d-f), we have A = 00 so that the NM leads do not play any role. In 
the first column, the length L\ of the first SC wire varies from 3 . 57 T to 7 . 5 - 7 T, while L 2 = 6 . 7 - 7 T and L 3 = 6 . 37 T are held fixed; 
Lj ~ 77 in all cases. In the second column, L\ varies from 14.5-7T to 18.5-7T, while L 2 = 12 . 3 - 7 T and L 3 = 12 . 77 r; hence Lj ~ 2 tj. In 
the third column, L\ varies from 70.57T to 74 . 57 T, while L 2 = 68 . 37 T and L 3 = 68 . 7 - 7 r; hence Lj 7$> 17 . The first row shows top 
views of surface plots of Gc as a function of L\ and E/A lying in the gap —A < E < A. The second row shows the energies 
E/ A of the sub-gap modes as a function of Li for a system with just three SC wires; the energies are found by the vanishing 
of the determinant of a matrix as explained in the text. The peaks in Gc shown in figures (a-c) match well with the sub-gap 
energies shown in figures (d-f). The third row shows surface plots of Gn 2 as a function of Li and E/A. 


In the regime of large wire lengths, the normal conduc¬ 
tance Gn 2 shows some interesting properties. From the 


discussion above, we would conclude that Gn 2 should be 
very small at this length regime. But in Fig. [2] (i) we 



















can see that Gn2 has discrete conductance peaks exactly 
at E = 0 at some particular values of L\, and Gjv 2 is 
quite high at those points. If we increase L\ more, the 
value of Gn2 decreases but the peaks still remain. This 
shows that the peaks are robust and are solely due to 
the three-wire geometry of the system as these kinds of 
peaks do not appear in a simple NM-SC-NM system^. 
We find that at the locations of these conductance peaks, 
the amplitudes t\j , t 2 j, and t&j take such values that 
t2je~ Lj ^, t^e Lj ^ and are signifi¬ 

cantly larger compared to their values when there is no 
conductance peak. Hence for these values of the ampli¬ 
tudes, the wave functions inside the SCs of the sub-gap 
modes are not negligible near the NM-SC junctions. So at 
these special points, an incoming electron from NM 1 can 
easily enter the SC and then get transmitted to the NMs 
on the other sides by coupling to these sub-gap modes. 

In the second column in Fig. [ 2 j we have taken L\ = 
14 . 57 T to 18 . 5 - 7 T, L2 = 12 . 3 tt and L3 = 12.77T. In this 
intermediate length regime, Li is of the order of the de¬ 
cay length £. Hence, the sub-gap modes will now hy¬ 
bridize with each other and their energy will split from 
zero. But the hybridization will be less than that in the 
Li < r) regime. Hence the energy splitting will also be 
less which is evident from the top view of Gc in Fig. [ 2 ] 
(b). In that plot we cannot see the six sub-gap modes 
separately. But there are states at non-zero energy due 
to the splitting, and the energy gap oscillates with L\. 
The highest value of Gc is lower than that in the Li 77 
regime as the incoming electron now has a finite proba¬ 
bility of transmitting to the other NM leads. 

The peaks in Gn2, as discussed in the large length 
regime, are present in this intermediate case also. But 
now Gn2 is non-zero at almost all values of L\ because of 
the appreciable hybridization among the sub-gap modes. 

Let us now compare our results with the energy plots 
we get for a SC box made of three SC wires. In these plots 
we use the same length regimes to compare it directly 
with our results. In Fig. [ 2 ](d), namely, when L\ < 77, we 
see that for a fixed L\ there are four points at different en¬ 
ergies where D = 0 . These four points correspond to the 
four sub-gap modes in Fig. [ 2 ] (a). We cannot see the other 
two modes because they are very close to E/A = ±1 and 
are therefore beyond our resolution. This difference be¬ 
tween the numerical results for the conductances and for 
the determinant occurs because in the conductance cal¬ 
culation the barriers at the NM-SC junctions have a finite 
strength while the determinant is calculated for a SC box 
with infinitely large barriers to the NM leads. 

As we move to the intermediate length regime, Fig. [2] 
(e) shows that the D = 0 lines approach zero energy; this 
occurs because the coupling between the sub-gap modes 
inside the SC box and therefore the splitting decreases. 
This is exactly what we see in the top view of Gc in 

Fig- HI (b). 

In the large length regime all the D = 0 lines stay 
almost exactly at E = 0 for all values of L\ as we can see 
from Fig. [ 2 ] (f). Now the sub-gap modes in the SC box 


are almost decoupled. This matches with the top view of 
Gc for large length, i.e., Fig. [ 2 ] (c). We conclude that our 
numerical results for Gc match well with the results we 
get for a SC box with infinite barriers at the three ends. 

3.2 One of the Aj’s with a different sign 

In this section we consider the case Ai <0 and 
A2, A3 > 0, but we take the magnitudes of all the Aj’s 
to be the same. As the electron is coming from the NM 1 
side, the results for Ai < 0 are different from the cases 
with A 2 < 0 or A3 < 0 . We will again calculate the 
conductances for three different length regimes and the 
determinant for a SC box. We then analyze the numeri¬ 
cal results to see the differences in the conductances com¬ 
pared to the case when all the Aj’s have the same sign. 
We will take the same values of the parameters as before, 
so that the length scale 77 = 20. 

In the first column of Fig.[ 3 l we take L\ = 2.5-7T to 6.5-7T, 
L2 = 5.77r and L3 = 5.37T. In this regime, from the top 
view of Gc in Fig. [ 3 ] (a), we can clearly see that there are 
four sub-gap modes in the system. The sub-gap modes 
have non-zero energies as the decay length £ is greater 
than the length of the SCs. Instead of three sub-gap 
modes at the junction of three wires, there is now only 
one sub-gap mode. In this length regime both Gc and 
Gn2 are finite as we can see from the plots. We have 
already discussed the reason behind this earlier when all 
the Aj’s have the same sign. 

In the third column, we choose L 1 = 70 . 57 T to 74.57T, 
L2 = 68.37T and L3 = 68.77 t; hence L\ 77 and is there¬ 
fore much larger than the decay length. Similar to the 
earlier case when all the A j ’s are positive, here we again 
get Gc almost equal to 2 (in units of e 2 /h) and Gn2 
close to zero as the sub-gap modes are almost decoupled. 
The energy splitting of the sub-gap modes is also close 
to zero. Instead of having different energies, now all the 
four sub-gap modes lie close to zero energy as we can see 
from Fig. [ 3 ] (c). Near the ends of the gap, i.e., E/A = ± 1 , 
Gn2 rises which is due to the single particle density of 
states which is non-zero at those energies. 

Unlike the case of all Aj’s positive, Gn2 now does not 
show any conductance peaks at E = 0 at any particular 
values of L\. We have looked at the numerical values of 
t±j, i 2 j, and t^j at E = 0 and at the values of L\ 
where Gn2 had peaks for the case of all Aj’s positive. 
We find that t 2 je~ Lj ^, and 

are quite large in this case also. However t n 3 and t a 3 
are almost equal; hence Gn 2 = \t n 3\ 2 ~ | t a 3 \ 2 is almost 
equal to zero. This is an interesting difference between 
the cases of all Aj’s having the same sign versus one of 
them having a different sign. 

For the intermediate regime we have chosen L\ = 14 . 57 T 
to 18 . 57 T, L2 = 12.77r and L3 = 12.77r, so that L\ ~ 77. In 
Fig. El (b) we see the that the sub-gap modes lie very close 
to E = 0 . There is a small energy splitting; this splitting 
is due to the finite coupling of the sub-gap modes as the 
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decay length is of the order of L\. We find that Gn2 is 
almost zero and very small compared to Gc- 

Next we consider a SC box and plot the determinant 
D of the matrix of amplitudes as a function of E/A and 
L\ as we did earlier, taking Ai < 0, and A2, A3 > 
0 . From the second row of Fig. [ 3 l we can see that the 


D = 0 lines match well with our numerical results for the 
conductance peaks in the different regimes of length. 

We have checked that instead of taking Ai < 0 , if we 
choose A 2 < 0 or A3 < 0 all the numerical results are 
qualitatively similar. 



FIG. 3 : Plots of conductances and energies of the sub-gap modes when Ai < 0 , while A2, A3 > 0 : the parameters chosen 
are kF = 1 , m = 0.5 and —Ai = A2 = A3 = A = 0 . 1 , so that rj = 20 . In figures (a-c) and (g-i), the strength of the NM-SC 
barriers is A = 5 , while in figures (d-f), we have A = 00 so that the NM leads do not play any role. In the first column, the 
length L 1 of the first SC wire varies from 3.57T to 7.57T, while L 2 = 6.7-7T and L3 = 6.3-7T are held fixed; Lj ~ r/ in all cases. In 
the second column, L\ varies from 14 . 5 - 7 T to 18 . 5 - 7 T, while L2 = 12.3-7T and L3 = 12.77T; hence Lj ~ 2 r/. In the third column, Li 
varies from 70 . 57 T to 74.57T, while L2 = 68.37T and L3 = 68.7-7T; hence Lj rf. The first row shows top views of surface plots 
of Gc as a function of Li and E/A lying in the gap —A < E < A. The second row shows the energies E/ A of the sub-gap 
modes as a function of L\ for a system with just three SC wires; the energies are found by the vanishing of the determinant 
of a matrix as explained in the text. The peaks in Gc shown in figures (a-c) match well with the sub-gap energies shown in 
figures (d-f). The third row shows surface plots of Gn2 as a function of L 1 and E/A. 
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We would like to point out that our results differ in 
many significant ways from those in earlier papers (such 
as Ref. I55T) which considered a single SC wire lying be¬ 
tween two NM leads. In a single SC wire, there are only 
two sub-gap modes lying at the ends of the wire; they 
give rise to two peaks in the conductance. If the p- wave 
pairing A changes sign somewhere inside the wire, two 
more sub-gap modes appear there; however these do not 
have an appreciable effect on the conductances. In con¬ 
trast to this, our system has several sub-gap modes, three 
at the NM-SC junctions and one or three at the junction 
of three SC wires depending on the relative signs of the 
Aj’s. All these sub-gap modes give rise to peaks in the 
conductance. To the best of our knowledge, this is the 
first model to be studied in which there are such a large 
number of conductance peaks whose locations precisely 
match the energies of all the sub-gap modes. 

4 Analytical understanding of Majorana 
modes at a junction of three SC wires 

In this section we will study how many zero energy 
Majorana modes appear near a junction of three SC wires 
if the SC wires are semi-infinite and the NM leads are 
absent. Namely, the system only consists of three semi¬ 
infinite SC wires which meet at Xj = 0. Therefore out of 
the four scattering amplitudes ty, t 2 j.tsj and ty, only 
those two will be non-zero for which the SC wave function 
V’scj —> 0 as Xj — > oo. From Eqs. (ED and ED, we know 
that ipscj is normalizable when only ty and t 2 j are non¬ 
zero. Keeping this in mind we can write down the wave 
functions for the semi-infinite wires. 

Ac, = til e"* 11 ’ ( A ) 

+ ( _ s 9 „(y) e -* ) ■ ( 2 3 ) 

We want to find Majorana modes which are at zero 
energy and are localized near the junction. For E = 0, 
e 1 ^ = [E — iy/ A 2 — E 2 )/A reduces to e 1 ^ = —i. Substi¬ 
tuting this in Eq. C3l) . we get 

’P.ci = as „(A,) ) 

( -, sgn(Aj) ) • ( 24 > 

Now we use the boundary conditions in Eqs. (1131) to find 


the Majorana modes. Using «m for the Cj 's we get 
iki tij — ik 2 t 2 j 

= Mji (tn + t 2 i) + Mj 2 (ti 2 +t 22 ) 

+ (ti 3 + <23) 

- (25) 

while for the dj’s we get, 

sgn(Aj) [iki ty — ik 2 t 2 j] 

= Mji sgn{Ai) (tn+t 2 i) 

+ Mj 2 sgn( A 2 ) (t 12 +t 2 2) 

+ Mj 3 sgn( A 3 ) (ti 3 + t 2 3 ) 

^AA (tij+t2j) (26) 

So we have a total of six equations, three for the Cj’s 
and three for the dj’s. From Eqs. E5l) and (14B1) it can be 
shown that if all the Aj’s have the same sign, the equa¬ 
tions for the c’s and d’s are exactly the same. So instead 
of six equations we get only three independent equations. 
But we have six variables in this problem since we have 
three wires and each of the wires has two different scat¬ 
tering amplitudes corresponding to k\ and k 2 . Hence we 
have three equations for six variables. Hence, there are 
three independent solutions and therefore three Majo¬ 
rana modes at the junction. 

Now we consider the case when one of the Aj’s has 
a different sign from the other two; let us assume that 
A 2 and A 3 have the same sign while Ai has the oppo¬ 
site sign. In this case, it is straightforward to show that 
tn + t 2 i = 0. So we now have five independent variables 
instead of six. Using the above relation it can be shown 
that the equations for c 2 and C 3 are the same as the equa¬ 
tions for d 2 and respectively. But the equations for Ci 
and di are different. So we have four independent equa¬ 
tions and five variables which implies that we have only 
one independent solution and hence only one Majorana 
mode. 

These are the only two combinations of the signs of 
Aj’s which will give different results. Eq. (1141) implies 
that any other combination will be similar to either the 
first case or the second case. To summarize, the two cases 
are as follows. 

(i) If all three Aj’s have the same sign then there are 
three Majorana modes at the junction. 

(ii) If one of the Aj’s has the opposite sign to the other 
two then there is only one Majorana mode at the junc¬ 
tion. 

The existence of one or three zero energy Majorana 
modes can also be understood as follows. A well known 
lattice model of a p -wave superconducting wire is the 
Kitaev chaini We first consider a single Kitaev chain. 
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The electron operators c n and c' n at site n can be writ¬ 
ten in terms of two Hermitian operators, called a n and 
bm as c n — (l/2)(<Xn T and c— (l/2)((2y 1 ib n ^j. 
(These operators satisfy the anticommutation relations 
{a m ,a n } = {b m ,b n } = 2 S mn ). The Hamiltonian of the 
Kitaev chain only has terms of the form ia m bm not ia m a,n 
or ibmbn ■ Such a Hamiltonian has an “effective time re¬ 
versal symmetry”, namely, it is invariant under complex 
conjugation of all complex numbers along with a n —> a n 
and b n —> — 6,^1. This symmetry implies that if we look 
at eigenstates with zero energy, their wave functions in¬ 
volve only the a n or only the b n , not both. A long Kitaev 
chain has zero energy Majorana modes at the two ends. 
Depending on the sign of the p -wave pairing, positive or 
negative, it turns out that the Majorana mode at the 
left end is of type a and at the right end is of type b or 
vice versasi. We now consider three semi-infinite Kitaev 
chains which meet at one site; the three Kitaev chains 
could have the same or different signs of the Aj’s. If all 
the Aj’s have the same sign, the three Majorana modes 
at the junction will all be of the same type, say, a. Then 
they will not mix with each other since the Hamiltonian 
has no terms of the form ia m a n ; hence they will remain 
at zero energy. On the other hand, if one of Aj’s has a 
different sign from the other two, two of the Majorana 
modes will be of one type, say a , and the other will be 
of type b. Now the b mode can mix with the two a’s; as 
a result the energies of two of the Majorana modes will 
become non-zero, namely EE, while one Majorana mode 
will stay at zero energy. (Written in terms of the a and b 
operators, the Hamiltonian of the system can be seen to 
have an £ —> —E symmetry; hence states can mix and 
move away from zero energy only in pairs). 

Having seen that a Kitaev chain has two kinds of Majo¬ 
rana end modes, called a and b, we can write down an ef¬ 
fective Hamiltonian which describes the sub-gap physics 
of our three-wire system. This Hamiltonian will only 
contain the operators for the Majorana end modes (the 
operators in the rest of the SC wires are not required) 
and the electron operators in the NM wires^ i 16 i 25 i 70 . The 
Hamiltonian will have three kinds of terms. First, the 
Hamiltonian at the junction of the three SC wires can 
be written as follows. If all three Majorana modes are 
of the same type, then no coupling between them is al¬ 
lowed due to the effective time reversal symmetry. But 
if two of the modes are of one type (say, type a on wires 
1 and 2) and the third mode is of the other type (b on 
wire 3), then this part of the Hamiltonian will take the 
form H i = i(w\ai + w 2 a 2 )b 3 , where w\ and w 2 are two 
couplings (assumed to be real). Next, on each wire j, 
the Majorana mode near the junction with the NM wire 
is of the opposite type to the mode near the junction 
of the three SC wires; namely, at the junctions with 
NM wires, the modes will be given by b \, b 2 and a 3 . 
The hybridization between the two modes at the ends of 
each wire will give rise to three more couplings given by 
H 2 = i(w 3 aibi + wia 2 b 2 + w 5 a 3 b 3 ). (The couplings w 3 , 
uq and w 3 go to zero exponentially as the lengths of the 


SC wires become large). Finally, we have couplings at 
the junctions of the SC and NM wires between the Ma¬ 
jorana modes at Xj — Lj — e and the electron operators 
in the NM wires at Xj = Lj + e (following the notation in 
Eq. ©)■ We can think of this coupling as arising from 
a hopping between the operator c sc = (1/2) (a + ib) at 
Xj = Lj — e in the SC wire and the operator c nm at 
Xj = Lj + e in the NM wire, namely, a term of the form 
E sc c nm + 4 m c sc . If the Majorana mode is of type a, we 
set b = 0 and get c\ c c nm + 4 m c sc = (l/2)a(c nm - c+ m ), 
while if the Majorana is of type b, we set a = 0 and get 
c\ c Cnm + cl m C sc = -( i/2)b(c nm + C^m)- We thuS § et a 
coupling of the form 

H 3 = iw 6 h [ c nm (Li + e) + 4 m (Iu + e)] 

+ iw 7 b 2 [c nm (L 2 + e) + 4 m (L 2 + e)] 

H” ^8 [077,777, (L 3 + e)-cL(L 3 + e)]. ( 27 ) 

The complete effective Hamiltonian is given by the sum 
of Hi, H 2 and H 3 and has eight parameters Wi. 

The parameters w , can be determined as follows. We 
first use the microscopic Hamiltonian in Eq. ©, the 6 - 
function potentials A at the NM-SC junctions and the 
matrix M at the junction of three SCs to compute the 
energies and wave functions of all the sub-gap modes. We 
then fit the values of these quantities with those calcu¬ 
lated from the effective Hamiltonian described above in 
order to find the values of w j. We will not carry out this 
exercise here. The effective Hamiltonian can then be used 
to compute the various conductances of the system-^. 

5 Renormalization group equations for a 
junction of several wires 

In this section, we will consider the effect of electron- 
electron interactions on the conductances of a system of 
two or more NM wires which meet at a junction where 
there is a SC region (we can consider this region to be 
a SC dot). This system differs from the one studied in 
the earlier sections in two major ways. First, the region 
of three SC wires considered in Sects. 2-4 will now be 
taken to be a single region. Further, the only role played 
by this region will be to give rise to a scattering matrix S 
for the NM wires. Due to the superconductivity, we have 
to consider the cases of both incident electrons and inci¬ 
dent holes, and both normal and Andreev reflection and 
transmission processes; S will therefore be a 2TV x 2 N ma¬ 
trix for the case of N NM wires meeting at the junction 
with the SC region. Second, we will only study the con¬ 
ductances at energies much larger than A, in contrast to 
the earlier sections where we only looked at the sub-gap 
conductances. The reason for considering energies much 
larger than A (but much smaller than the band width of 
the system) is that the analysis given below only works 
when the conductances are slowly varying functions of 
the energy; this will become clearer as we proceed. We 
will take the interactions to be present only in the NM 
wires and derive RG equations which will tell us how the 
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S'-matrix evolves when we start at the length scale of the 
SC region and then increase the length scale to expand 
from that region into the NM wires. 

In one dimension, the technique of bosonization is 
often used to study systems with short-range density- 
density interactions^; in this method a system of inter¬ 
acting fermions is mapped to a system of non-interacting 
bosons. This method has the advantage that interac¬ 
tions of any strength can be dealt with. However this 
method runs into difficulties in the presence of junctions 
of three or more wires and superconductivity for the fol¬ 
lowing reasons. As we have seen, a junction is charac¬ 
terized by a matrix M which relates the electron fields 
on different wires in a linear way. Since bosonization 
relates electron operators to exponentials of bosonic op¬ 
erators, a linear relation between fermionic fields in dif¬ 
ferent wires translates, in general, to a non-linear relation 
between the bosonic fields; this makes it difficult to use 
bosonization. (However, there are some special forms of 
the junction matrix when bosonization can be applied. 
These special cases correspond to the magnitudes of all 
the reflection or transmission amplitudes being equal to 
either zero or 1, namely, either perfect reflection and 
no transmission, or no reflection and perfect transmis¬ 
sion^). Next, bosonization works best if the system is 
gapless and the energy-momentum dispersion is linear for 
both the fermionic and the bosonic theories which are re¬ 
lated to each other. However superconductivity produces 
a gap proportional to the SC pairing A, and the disper¬ 
sion is not linear for energies of the same order as A. So 
bosonization works only if we treat A as a perturbation 
(as was done in Ref.114 for example). 

We will therefore use a different approach which di¬ 
rectly uses the fermionic language^—. Unlike bosoniza¬ 
tion, this approach is useful only if the interaction 
strengths are weak in all the wires. However it has the 
advantage that it works for any form of the scattering 
matrix which characterizes the junction. The results ob¬ 
tained by this method and those obtained by bosoniza¬ 
tion will of course match in the parameter regimes where 
both methods work, namely, when the interactions are 
weak and the junction scattering matrix has some spe¬ 
cial forms. We will now describe this method in detail. 

We begin with a second quantized fermionic field c(x) 
and the corresponding hole field d(x), where d(x) = 
c' ( x ), for a single semi-infinite NM wire which goes from 
x = 0 to oo and is connected to a SC region at x = 0. 
At low temperatures, only low-energy processes are of in¬ 
terest and these only involve modes near the Fermi mo¬ 
menta ±fcj f; we will therefore consider only these modes. 
We therefore write the second-quantized field c(x) as 

c{x) = a{x)e- ikFX + c 0 (x)e ikFX , (28) 

where cj and Co denote the fields of incoming and out¬ 
going electrons respectively. We take these fields to be 
slowly varying on the length scale of 1 /kp as we have sep¬ 
arated out the rapidly varying factors e ±lkpx . Namely, 
the fields cj, cq have momentum components k such that 


| A: | -C &F- We can then use a linear approximation for the 
dispersion relations of these fields so that E = ±vpk for 
ci and co respectively, with vf being the Fermi velocity. 

Similarly, for the second-quantized field d(x) we write 

d(x) = d!{x)e ikFX + d 0 (x)e~ lkFX . (29) 

Note that the rapidly varying exponential terms multi¬ 
plying di and do are the opposite of those multiplying cj 
and co ■ This is because destroying an electron is equiva¬ 
lent to creating a hole, so that Eh = —E e where Eh and 
E e are the energies of the electron and hole respectively. 

Let us now introduce a short-ranged density-density 
interaction between the electrons of the form 

Hint = \ J J dxd V p(x)V(x - y)p(y). (30) 

We assume that V(x) is a real and even function of x. 
The density p is a function of the second quantized fields 
given by p(x) = c'(x)c(x) = d{x)d'{x) = —d'{x)d{x) 
(using the anticommutation property of the fermionic 
fields). Using Eq. (f28l) . we obtain the expectation val¬ 
ues 

< p(x) > = < c|c/ > + < CqCo > 

i ^ t „ \ 2 ikpx , t _ ^ ^—2 ik&x 

+ < CjCo > e + < c' 0 ci > e . 

(31) 

Next, we will assume that V(x) is so short-ranged 
that x and y, which appear as arguments of the den¬ 
sity fields, can be set equal to each other except when 
the corresponding term in Hi nt becomes zero. Using this 
assumption and the anticommutation relation between 
the fermionic fields, we get 

Hint = g J dx 4cr4,co, (32) 

where g is related to the Fourier transform of V(x) as 
g = V’(O) — U(2A:f). From this expression it is clear 
that g = 0 if V(x) is a ^-function. Hence V(x) must 
have a finite range for the interaction to have an effect. 
For each NM wire, the interaction is described by the 
single parameter g. The value of this parameter may be 
different for different wires which we will denote by gj on 
wire j. We define a dimensionless quantity ay as 


where we assume vf to be the same in all the NM wires. 

Let us briefly discuss how the interaction parameter 
appears in the formalism of bosonization. For spinless 
fermions, which is relevant here as we are studying p -wave 
SCs, the bosonic theory is characterized by two quanti¬ 
ties, the velocity v of the excitations and a dimensionless 
parameter K (called the Luttinger parameter) which is a 
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measure of the strength of the interactions between the 
fermions. These are related to v F and a as 


= v F (l - a 2 ) 1/2 , 
1/2 

K = 


1 + a 


(34) 


Note that K = 1 if a = 0 (non-interacting fermions), 
while K < 1 ( K > 1) for a > 0 (a < 0), namely, repulsive 
(attractive) interactions. For weak interactions we get 
v = Vp and K = 1 — a to first order in a. We will do our 
RG analysis in the limit that a is small and positive in 
each wire. 

As we have discussed earlier, it is generally difficult 
to bosonize a system with junctions. We will therefore 
use a different method which will enable us to derive 
RG equations directly for the scattering matrix of the 
junction. As we will see, this method only works up to 
first order in the interaction parameters ay. The basic 
idea of this method is the following. In the presence of 
non-zero reflection amplitudes rjj at the junction, the 
density of non-interacting fermions in the NM wire j will 
have Friedel oscillations with wave number 2 k F . When 
an interaction is turned on, an electron can scatter to an 
electron or a hole (Andreev reflection) from these oscil¬ 
lations with an amplitude which is proportional to the 
parameter aj . Ref. |8J used this idea to derive RG equa¬ 
tions for an arbitrary S-matrix describing the junction of 
two semi-infinite wires. An RG analysis was then done 
for junctions of more than two wires, without supercon¬ 
ductivity in Refs. [85| and [8(1 a nd with a s-wave super¬ 
conducting junction in Ref. [83- We will carry out an RG 
analysis for our system where several NM wires meet at a 
junction with a p- wave SC region. We expect the results 
to be much richer than those for a junction of NM wires 
when no superconductivity is present. 

We will begin our analysis by deriving the form of the 
density oscillations in one particular NM wire close to the 
junction with the SC region. We will consider separately 
the cases of an electron and a hole coming in from a NM 
lead. 


5.1 Processes related to an incoming electron 


of k in Eq. (13511 . Although we are only interested in values 
of k close to zero, it is mathematically convenient to take 
the range of k to be — oo to oo even though the range is 
finite and given by the band width in real systems. The 
expectation value of p[x) in terms of c k (x ) is given by 

rib 

< p{x) > = J 7 ^z c* k c k , (36) 

where we have taken the lower limit of k to be — oo. We 
see that < p > has a constant piece po which can be 
eliminated by normal ordering. We are then left with 


< P {x) > — po = ^< e e~ 2 ' kFX - r ee e 2ikpx ). (37) 

It is clear that these terms arise entirely due to the 
interference between the incoming and reflected waves. 
Substituting Eq. (1371) in Eq. (1311) we see that < c\ci + 
CqCq >= Po is a constant, while 


t ir ee 

< C ° C1 > = w 

t t * 

< c}co > = < CqCi > 


ir ee 

4-7TX 


(38) 


An important point to note here is that there will also 
be a contribution to c(x) and therefore to < p{x) > from 
the waves which are transmitted from the other wires. 
As a wave transmitted to one wire from any other wire 
is incoherent with the incident and reflected waves of 
the first wire (we are assuming that waves incident from 
different NM leads are phase incoherent with respect to 
each other), there is no interference between these waves. 
Since the waves transmitted from the other wires only 
contribute to an outgoing wave in this wire, there is no 
interference and their contribution to < p{x) > is inde¬ 
pendent of x. Hence, it can be absorbed in po . We con¬ 
clude that the Friedel oscillations in Eq. (1371) in a given 
wire arises only from the reflections within that wire. 

Now, in our system we have both normal and Andreev 
reflections. So there will be some non-zero expectation 
values for the operators which connect electrons and holes 
such as c\do and d^cj. Using Eq. (1551) . we find that the 
expectation values of these operators are given by 


An incoming electron can be either 

(i) normally reflected with amplitude r ee , or 

(ii) Andreev reflected to a hole with amplitude rh e . 

For momenta near k F we can write the wave functions 
for electrons and holes as 


< c\do > 

< dpCi > 


J —OO 

• 3 k 

lr he 
47 TX 


dk 

2 ir 


Thee- 


i2kx 


irhe 
47TX ’ 


(39) 


c k {x) = Cl e~ lkFX + c 0 e ikFX 

— e ~i(k F +k)x i(k F +k)x 

T ' ee c 5 

d k {x) = doe~ lkFX = r he e-^ kF ~ k ^, (35) 

where |fc| <C k F . In the ground state of a non-interacting 
system, all the energy states below the Fermi energy E F 
are filled for electrons; this corresponds to negative values 


To evaluate the integral in the first equation in (l39l) . we 
must introduce a factor like e tk which cuts off the contri¬ 
bution from the lower limit k —> —oo, and we then take 
the limit e —> 0. Further, we have assumed that rh e varies 
slowly with k so that it is a reasonable approximation to 
take it outside the integral over k in the first equation. 
This is the reason why our analysis only works at ener¬ 
gies which lie far from the SC gap; for those energies the 
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reflection and transmission amplitudes are slowly varying 
functions of the energy. In contrast to this, the sub-gap 
conductances have sharp peaks due to various sub-gap 
modes; hence the reflection and transmission amplitudes 
are not slowly varying functions of the energy if it lies in¬ 
side the SC gap. We note that a renormalization group 
study at energies within or close to the SC gap, where 
the reflection and transmission amplitudes vary rapidly, 
has been carried out in Ref. |90|. 

Next, we derive the reflections of the electrons and 
holes from the Friedel oscillation by using the Hartree- 
Fock decomposition of the Hamiltonian in Eq. (1321) . We 
have 

n OO 

Hint =-g dx [< c\co > CqCi + < CqCi > c\co 

Jo 

+ < c\cq > cicq + < cicq > c\cq], (40) 


Using the expectation values of the various operators 
derived earlier and the identities c/co = c/dj), = — d' 0 Ci 
and c\cq = c\do , we obtain 

gf°°dx t . * t 

H mt = — / — [ir ee c' 0 ci - ir ee c\c 0 

4tt Jo x 

- ir he d[ ,c/ + ir* he c\d 0 \■ (41) 

We can now derive the amplitude to go from an in¬ 
coming wave to an outgoing wave under the action of 
exp (—iHintt). We begin with an incoming electron with 
momentum k. Various processes can now occur. 

(i) An incoming electron with momentum k can go to an 
outgoing electron with momentum k’ under the action of 
Hi n t ■ The corresponding amplitude is 


/ lik 1 

2 nS(E k - E k ,) 

x < outgoing,e,k'\Hi n t\incoming,e,k > 

f°° — r P~ 2ikx 

^V F J 0 X 


(42) 


To obtain the above expression we have used Eq. (HOI) , the 
dispersion relation E = v F k (which implies 5(E k —E k i) = 
(1 /v F )d(k—k')), and the wave functions e ±l ( fcF+fc ) x of the 
outgoing and incoming electrons respectively. 

To derive RG equations for quantities like r ee from 
Eq. (l42l) . we will integrate x over a small interval going 
from L to L + dL = Le dl . Here l is the logarithm of the 
length scale, and we can write l = ln(L/a ), where a is 
a short distance scale (which is the size of the supercon¬ 
ducting region forming the junction) from which we will 
begin to integrate the RG equations. Eq. (l42l) then gives 


dl = ^ dl , (43) 

Attv f 2 

where we have used Eq. (l33l) . 

(ii) Similarly we find that the amplitude to go from an 
outgoing electron to an incoming electron is given by 


ar 


* 

ee 


di. 


(44) 


(iii) Due to the presence of the SC region, Andreev reflec¬ 
tion can also occur, namely, an incoming electron can go 
to an outgoing hole under the action of Hi nt as given in 
Eq. (fllip . We can calculate the amplitude of this pro¬ 
cess in the same way as we did for normal reflection 
above. The amplitude to go from an incoming electron 
with momentum k to an outgoing hole with momentum 
k' is found to be 


ar he 

2 


dl. 


(45) 


(iv) The amplitude to go from an outgoing hole to an 
incoming electron is given by 


rvv* 

-J^- dl. (46) 

This completes the list of processes which can occur if 
we start with an incoming electron. 


5.2 Processes related to an incoming hole 


All the processes we discussed in the previous section 
can be studied if we start with an incoming hole. An 
incoming hole can be either 

(i) normally reflected to another hole with amplitude r^, 
or 

(ii) Andreev reflected to an electron with amplitude r e h- 
For momenta near k F we can write the hole and elec¬ 
tron wave functions as 


d k {x) = d ie lkFX + d 0 e~ lkFX 

— e i(k F -k)x _|_ rhhe -i(k F -k)x ^ 

c k (x) = c 0 e lkFX = r eh e^ kF+k ^. (47) 

For holes, we can use Eq. (1M1) and the fact that cU = 
ddJ = —dJd (plus a constant) to write the ground state 
expectation value of p(x) as 

r° rjb 

< P(x) > = - ^^44 (48) 


We can then find expectation values of various operators 
following a procedure similar to the previous section. We 
obtain 


< d\do > 

< dpdi > 

< Cpdi > 

< d\co > 


ifhh 

Anx 


' hh 


47 TX ’ 

* A 

lr eh 

Anx ’ 
ir e h 
4nx 


(49) 


Using these expectation values we find the amplitudes of 
various processes. 

(i) Under the action of Hi nt , an incoming hole goes to an 
outgoing hole with an amplitude 


otr h h 

2 


dl. 


(50) 


2 
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(ii) An outgoing hole goes to an incoming hole with an 
amplitude 


- dl. (51) 

(iii) An incoming hole goes to an outgoing electron with 
an amplitude 


ar eh 

2 


dl. 


(52) 


(iv) An outgoing electron goes to an incoming hole with 
an amplitude 


ar 


* 

eh 


2 


dl. 


(53) 


5.3 RG equations for the reflection and transmis¬ 
sion amplitudes 

The amplitudes of all the processes that we derived 
earlier can now be combined along with the S'-matrix at a 
junction of several NM wires with a SC region to calculate 
corrections to the S'-matrix. We will calculate all the 
corrections to first order in the interaction parameters 
ctj. 

We first consider corrections to the reflection ampli¬ 
tude r ee jj on wire j. To first order in otj, this gets 
contributions from the following processes. An incom¬ 
ing electron on wire j can 

(i) become an outgoing electron on the same wire by scat¬ 
tering from the Friedel oscillations in the density with the 
amplitude given in Eq. (1431) . 

(ii) reflect from the junction with amplitude r ee jj to be¬ 
come an outgoing electron, then become an incoming 
electron due to the Friedel oscillations according to (l44l) , 
and finally reflect from the junction as an electron with 
amplitude r ee jj. 

(iii) reflect from the junction with amplitude r ee jj to 
become an outgoing electron, become an incoming hole 
due to Friedel oscillations according to (1551) . and then 
Andreev reflect from the junction as an electron with 
amplitude r e h,jj- 

(iv) Andreev reflect from the junction with amplitude 
fhe.jj to become an outgoing hole, become an incoming 
electron due to Friedel oscillations according to (l46ll , and 
then reflect from the junction as an electron with ampli¬ 
tude r ee jj. 

(v) become an outgoing hole by Andreev reflection from 

the junction with amplitude rh e ,jj, become an incoming 
hole due to Friedel oscillations according to and 

then Andreev reflect from the junction as an electron 
with amplitude r e hjj- 

(vi) transmit through the junction to wire i (with i ^ j) 
as an electron with amplitude t ee ^j, turn from an out¬ 
going electron to an incoming electron on wire i due to 
Friedel oscillations according to (1441) . and then transmit 


through the junction to wire j as an electron with am¬ 
plitude t ee ji. 

(vii) transmit through the junction into wire i as an elec¬ 
tron with amplitude t ee ^j , turn from an outgoing electron 
to an incoming hole on wire i according to (1531) . and then 
transmit to wire j as an electron with amplitude t e h,ji- 
(viii) transmit through the junction into wire i as a hole 
with amplitude th e ,ij , turn from an outgoing hole to an 
incoming electron on wire j according to (|46D . and then 
transmit to wire i as an electron with amplitude t ee ,ji- 
(ix) transmit through the junction into wire i as a hole 
with amplitude th e ,ij , turn from an outgoing hole to an 
incoming hole on wire j according to ED, and then trans¬ 
mit to wire j as an electron with amplitude t e h,ji- 

Collecting all these terms, we find that the correction 
to r ee jj is given by 


dr, 


ee,jj 


d^eejj dl j 


Bee,jj ~ r,[ a j r ee,jj a j\ r ee,jj\ r ee,jj + Otj\r e h,jj \ T ee ,jj 
+<%j\ r he,ji | Tee,jj ~ a j r he,jj1~hh,jj r eh,jj 
+ 'y J {~ a ilee 1 ji r ee ^iitee,ij T a i^ee,ji r e h,ii^eh,ij 

+ a ithe,ji r he,iitee,ij ~~ a i^he,ji1~hh,ii^eh,ij)]- 

(54) 


Similarly, the transmission amplitude t ee ji from wire i 
to wire j can get corrections from the following processes. 
The incoming electron on wire i can 

(i) get reflected from the junction as an electron with am¬ 
plitude r ee: a, then become an incoming electron accord¬ 
ing to (1441) , and finally transmit to wire j as an electron 
with amplitude t eeji . 

(ii) reflect from the junction as an electron with ampli¬ 
tude r ee ,ii, become an incoming hole according to (1531) . 
and then transmit to wire j as an electron with amplitude 

teh,ji ■ 

(iii) Andreev reflect from the junction as a hole with am¬ 
plitude The,a, become an incoming electron according to 
(USD, and then transmit to wire j as an electron with 
amplitude t ee ji. 

(iv) Andreev reflect from the junction as a hole with am¬ 
plitude rhe,a, become an incoming hole according to (1511) . 
and then transmit to wire j as an electron with amplitude 

teh,ji ■ 

(v) transmit to wire j as an electron with amplitude t ee ji , 
become an incoming electron on wire j according to (1441) . 
and then reflect from the junction as an electron with 
amplitude r ee ,jj. 

(vi) transmit to wire j as an electron with amplitude 
tee,jii become an incoming hole on wire j according to 
ED, and then reflect from the junction as an electron 
with amplitude r e h,jj- 

(vii) transmit to wire j as a hole with amplitude th e ,ji , 
become an incoming electron on wire j according to (1461) , 
and then reflect from the junction as an electron with 
amplitude r ee ,jj. 
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(viii) transmit to wire j as a hole with amplitude th e ,ji , 
become an incoming hole on wire j according to (1511) . 
and then reflect from the junction as an electron with 
amplitude r e hjj- 

(ix) transmit to wire k (with k ^ i, j) as an electron with 
amplitude t ee ,ki, become an incoming electron according 
to (El) , and then transmit to wire j as an electron with 
amplitude t ee jk- 

(x) transmit to wire k as an electron with amplitude 
tee,ki> become an incoming hole according to (1551) . and 
then transmit to wire j with amplitude t e hjk- 

(xi) transmit to wire k as a hole with amplitude th e ,ki , 
become an incoming electron according to (1461) , and then 
transmit to wire j with amplitude t ee ,jk- 

(xii) transmit to wire k as a hole with amplitude th e ,ki , 
become an incoming hole according to m, and then 
transmit to wire j with amplitude t e hjk- 

Hence the total correction to t ee ji is given by 


dt 


ee,ji 


Bee,ji dl) 




ee,ji 


2 1 |^ee,£i | H” &iteh,ji r f'ee,iiTeh,ii 

~\~CYit ee ^ji \The 

~~ a j\ r ee,jj\ teeji + &j\r e h,jj\ tee,ji 

~kOijr ee jjr he jjth e ji — OtjVehJjT^jjtheJi 

+ E ( &ktee ,jk^ ee 

~\~Olkt e h,jk 1 " e h,kk^ee,ki + Oiktee,jk r he,kk^he,ki 

~ a -kteh,jk1'hh,kk^he,ki)]- (55) 


Similarly, we can find the corrections to 
all the other entries of the S'-matrix, namely, 

khli.jj ■ k he.jj : T eh, jj 5 tflh.j'/. • t/ie.ji and t'eh.ji ■ We 

now consider all these components of S to be functions 
of a length scale L , where L can vary all the way from 
a short distance scale a to a large distance scale which 
may be either the length of the NM wires or the thermal 
coherence length as discussed in Sect. 6.2. L and l are 
related as l = ln(L/a ), so that l = 0 when L = a. 
Eos. (15411551) then give us the RG equations 


d r ee,jj _ D 

dl ~ ™ 

dteeji r-j 


(56) 


Eqs. (1561) can be written in a more compact way. Given 
the matrix S and the parameters ctj , we define a matrix 
F whose non-zero elements are 


F — 

1 ee,jj ~ 

2 a j r ee,jji 


Feh,jj = 

1 

~ 2 a j re h,jji 


Fhe,jj = 

1 

~ 2 a j rhe ,jj’ 


Fhh,jj = 

1 

2 a j r hh,jj- 

(57) 


We can then show that the RG equations for the different 
elements of S (such as r ee jj and t ee ji given in Eq. (1561) 4 
can be written compactly in the form of a matrix equa¬ 
tion 


— = F-SF^S. (58) 

This is the central result of this section. Note that these 
equations are first order in the interaction parameters ctj 
since we have only considered processes with one scat¬ 
tering from the Friedel oscillations; hence the equations 
are valid only for weak interactions. (The interaction pa¬ 
rameters ctj do not themselves flow under RG since we 
are considering a spinless system^I). 

We can verify from Eq. (1551) that S continues to remain 
unitary under the RG flows; it also remains symmetric if 
it begins with a symmetric form since F is always sym¬ 
metric. We also note that the form of Eq. (1551) remains 
unchanged if S is multiplied either from the left or from 
the right by a diagonal unitary matrix with entries of the 
form 


Un=e^, (59) 

where the real numbers (f>j are independent of the length 
parameter l. Hence the fixed points discussed below 
will also remain unchanged under such phase transforma¬ 
tions. We will not distinguish between S'-matrices which 
differ only by such phase transformations since exper¬ 
imentally we generally measure scattering probabilities 
rather than amplitudes. 

6 RG fixed points, stability analysis and 
conductance 

6.1 RG fixed points 

We will now study the RG flows and find the fixed 
points and their stabilities. For a system with a junction 
of N wires, S' is a 2N x 2 N matrix which relates the 2 N- 
dimensional incoming and outgoing fields as (co, do) T = 
S(cj, di) T . S satisfies the following properties. 

(i) Unitarity: conservation of the particle current implies 
that S^S = I 2 N where I 2 N is the 2 N x 2N identity ma¬ 
trix. 

(ii) Time reversal symmetry: under time reversal, we 
complex conjugate all numbers, change t — > —t, c —»■ c*, 
and d —> —d*. From Eqs. (l35l) and (1471) . we see that this 
transforms cj —> Cq, co —> Cj, di —> —dp, do —► — d}. 
This implies that S satisfies S* = r 2 S*r 2 , where 

t 2 = ^ Iq j ^ , and Ijy is the N x N identity matrix. 

We therefore get S T = t z St z . 

If we want to find the fixed points analytically, we can 
use Eq. (1551) . dS/dl = F — SF^S = 0, which implies that 

SF* = F&. 


(60) 
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We can use this condition along with the two properties 
of S mentioned above. After finding the fixed points of 
the RG equations, we can study their stabilities. To do 
this, we write a fixed point of the S'-matrix as So, and a 
small deviation from this as eSi, where e is a small real 
parameter and Si is a matrix, namely, 

S = S 0 + eSi. (61) 

For a given So, we can find the various flow ‘directions’ 
Si such that Eq. (1551) takes the form 


de 

dl 


= Pe, 


(62) 


where f3 is a real number. The solution of this equation is 
e(l) = exp(/3^)e(0) where e(0) is given by the deviation of 
S from S 0 at the short distance scale a. We see that /? < 0 
indicates that S is stable against a perturbation in the 
direction of the corresponding Si, while /3 > 0 indicates 
an instability in the direction of Si. The case /3 = 0 
describes a marginal direction; for instance, this arises if 
we perturb a fixed point by a phase transformation as in 
Eq. (1591) which maintains it as a fixed point. 

It turns out to be difficult to find all the fixed points 
and their stabilities analytically, particularly for the 
three-wire case. We will therefore study this problem nu¬ 
merically by starting with some randomly chosen matrix 
S, evolve it according to Eq. (1551) . and see where it flows. 
However, we would like to begin with a matrix which sat¬ 
isfies the conditions S = I2N (which implies non-linear 
constraints on the elements of S) and S T = t z St z . To 
generate such a matrix, let us write S as S = exp(iA), 
where A is also a 2N x 2N matrix. Then the unitarity 
of S implies = A, and S T = t z St z then implies that 
A T = t z At z . These provide linear relations between 
the elements of A and are therefore easier to implement. 
(One can show that N(2N + 1) real parameters are re¬ 
quired to specify such a matrix A). After randomly gen¬ 
erating a matrix A which satisfies these conditions, we 
take S = exp(iA) as the starting matrix for the RG flows 
and see where it flows at large distance scales. 

We first discuss a two-wire problem as a guide to the 
more complicated three-wire problem. Here we have two 
NM wires with a SC region present at their junction. 
Then S is a 4 x 4 matrix of the form 


/ r ee 

tee 

Teh 

teh 

tee 

Tee 

teh 

Teh 

The 

the 

Thh 

thh 

V the 

The 

thh 

Thh 


(63) 


(The condition S T = t z St z implies that The = —r e h and 
the = — t e h )■ Here we have assumed for simplicity that all 
the scattering amplitudes are the same for wires 1 and 2, 
for instance, r ee ,n = fee ,22 and so on; these amplitudes 
can be complex in general. We will also assume that 
the interactions parameters are equal on the two wires, 

ai = 02 - 




FIG. 4: Typical RG flows of various scattering probabilities 
for two wires meeting at a junction with a p-wave supercon¬ 
ducting region. Figure (a) shows the case with all ay = 0.2, 
while (b) shows the case with all ay = —0.2. 


We will first list the fixed points that we have found 
analytically for S'-matrices of the form in Eq. (l63l) . We 
will present the values of all the reflection and transmis¬ 
sion amplitudes, with the understanding that two fixed 
points in which the amplitudes only differ by some phases 
as discussed in Eq. (l59l) will not be considered as differ¬ 
ent fixed points. We have found the following fixed points 
which are similar to the ones found for a junction with a 
s-wave superconductor in Ref. [87l 

(i) r ee = rhh = 1, and all the other amplitudes are zero. 

(ii) r e h = —The = 1, and all the other amplitudes are 
zero. 

(iii) t ee = thh = 1, and all the other amplitudes are zero. 

(iv) t e h = — the = 1, and all the other amplitudes are 
zero. 

(v) r ee = thh = r he = the = 1/2, and t ee = r h h = r eh = 
t eh = -1/2. 

(vi) t ee — thh = cos 0, t e h = —the = sin 8 , and all the 
other amplitudes are zero. Here 9 can be any real number 
from 0 to 27 t. We thus have a continuous family of fixed 
points labeled by 9. This is in contrast to the fixed points 
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(i), (ii) and (v) given above which are all discrete fixed 
points. We note that the fixed points in (iii) and (iv) 
are special cases of (vi) corresponding to 9 = 0 and n/2 
respectively. 

A simple way of finding fixed points analytically is to 
set all the reflection amplitudes equal to zero, r ee ,jj = 
r hhjj = r eh jj = r he jj = 0. Then the matrix F 
in Eq. m vanishes, and Eq. (1551) straightaway gives 
dS/dl = 0. All the other amplitudes in S are then con¬ 
strained by unitarity. This is how the family of fixed 
points in (vi) was found. 

We will now numerically study the fixed points and 
their stabilities for Eq. (1551) . We begin by randomly gen¬ 
erating a matrix A satisfying the conditions A = A and 
A T = t z At z . Given an A, we construct S = exp(iA), 
and then use Eq. (1551) to evolve S. In Figs. 0] (a) and 
(b), we show the typical RG flows of a number of reflec¬ 
tion and transmission probabilities, such as R ee = \r ee j 2 , 
Rhh = \rhh\ 2 , etc. We find numerically that for a gen¬ 
eral starting point, S always flows to the same fixed 
point which is described below. (By a general starting 
point, we mean that we do not start precisely at a special 
set of matrices which constitute some other fixed points 
which are not completely stable). For the case where 
ai = a 2 = 0.2 (i.e., the interactions are repulsive), the 
plots in Fig. [4] (a) show that S flows to a fixed point 
where R ee = Rhh = 1, and all the other probabilities are 
zero. If we start the RG evolution exactly at this point, 
S does not flow at all and remains at that point. Hence 
this is the fixed point in the case of repulsive interactions. 
This is a completely stable fixed point because the RG 
flows always approach this fixed point independent of the 
S that we start from; further, if we deviate a little bit in 
any direction from this fixed point, the RG flows take us 
back to the fixed point. 

Interestingly, we find even if none of the symmetries 
are present, namely, au ^ a.i (but both are positive), 
there is no time reversal symmetry ( S T ^ t z St z \ and 
the starting 5-matrix is not symmetric between the two 
wires, S always flows to a completely stable fixed point 
which is symmetric. More carefully speaking, one finds 
at the fixed point that although the phases of r ee jj and 
r hh,jj are generally not equal to each other (the phases 
depend on the starting value of S and the values of ctj), 
all their magnitudes are equal to 1. We thus have a 
symmetry restoration at the fixed point. 

For the opposite case in which a.\ = 02 = —0.2 (i.e., 
attractive interactions), the plots in Fig.[]](b) show that 
5 generally flows to a different fixed point where R e h = 
Rhe = 1, and all the other probabilities are zero. This is 
the completely stable fixed point in this case as we always 
reach this point no matter where we start (unless again 
we start precisely at some special set of matrices which 
are fixed points which are not completely stable), and 
small deviations from this fixed point in any direction 
also flow to zero. Once again, we find that the phases 
of r e hjj and rh e ,jj are generally not equal to each other 
but all their magnitudes are equal to 1 at the fixed point, 


even if there is no symmetry in the starting value of S 
and of the ay (provided that they are both negative). 

We now observe that since Eq. (1551) is linear in the 
ctj, the equation depends only on the combination otjl. 
Hence the direction of the RG flows for ctj > 0 and l —> 00 
is just the opposite of the flows for ctj < 0 and l —> 00. 
We therefore conclude that the stable fixed points for at¬ 
tractive interactions are exactly the same as the unstable 
fixed points for repulsive interactions. So the fixed point 
with R e h = Rhe = 1 is the completely unstable fixed 
point for repulsive interactions, i.e., if we start slightly 
away from this point in any direction, the RG flows will 
always take us further away from the point. 

We note that this numerical way of finding fixed points 
cannot detect fixed points which are partially stable and 
partially unstable, namely, stable in some directions and 
unstable in other directions. If we start near such a fixed 
point, we will generally always flow away from it regard¬ 
less of whether we take all the ctj to be positive or nega¬ 
tive. 

Finally, let us discuss the RG fixed points and their 
stabilities for the three-wire problem. We again assume 
for simplicity that there is complete symmetry between 
the three wires, both in the 5-matrix and the interactions 
strengths. 5 is therefore a 6 x 6 matrix of the form 


/ 

Tee 

tee 

tee 

Teh 

teh 

teh 

\ 


tee 

Tee 

tee 

teh 

Teh 

teh 


tee 

tee 

r ee 

teh 

teh 

Teh 



The 

the 

the 

Thh 

thh 

thh 



the 

The 

the 

thh 

Thh 

thh 

/ 

V 

the 

the 

The 

thh 

thh 

Thh 


We now present the results we find numerically, by 
starting with a randomly chosen matrix S = exp(*A) 
which has all the desired symmetries and then evolving it 
using Eq. (1551) . The typical RG flows of the various scat¬ 
tering probabilities are very similar to the ones shown 
for a two-wire system in Figs. [5] (a) and (b) for all the 
ctj equal to 0.2 (repulsive) and —0.2 (attractive inter¬ 
actions) respectively. We see that the completely sta¬ 
ble fixed point for repulsive interactions is again given 
by R ee = Rhh = 1, and all the other probabilities are 
zero. As we discussed above, the stable fixed points 
for attractive interactions, i.e., ctj < 0, are the unsta¬ 
ble fixed points for repulsive interactions. Using this 
property we find that the completely stable fixed point 
for attractive interactions, and therefore the completely 
unstable fixed point for repulsive interactions, is given 
by R e h = Rhe = 1, and all the other probabilities are 
zero. All these statements remain true even if the ctj 
are not equal to each other (although they must all have 
the same sign) and even if the starting value of S has 
no symmetries. Thus the completely stable and com¬ 
pletely unstable fixed points are similar for the two-wire 
and three-wire systems. 

Before ending this section, we would like to mention 
the work done in Ref. [70| on a junction of a supercon¬ 
ducting wire and two non-superconducting wires where 
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there are interactions. A non-trivial stable fixed point 
was found there which has perfect normal reflection for 
one linear combination of the electron fields in the two 
wires and perfect Andreev reflection for the other linear 
combination. This fixed point occurs because Ref. Iz 3 
considers an interaction between the wires which mixes 
the two electron fields. We have not considered such 
inter-wire interactions in our model and therefore do not 
find such a non-trivial fixed point. 


the length scale where the RG flows must be stopped is 
L w . For a general perturbation around the stable fixed 
point, we have \r he \ 2 ~ e~ 4al , \t he \ 2 ~ e~ 2a \ and \t ee \ 2 ~ 
e _2 “ J , where l = ln(L w /a). We then find that Gc = 
2 (Ke | 2 + 2 |fhe| 2 ) scales as 

G c ~ Cl e- 4 “* + c 2 e~ 2al 

~ Cl + Cz ('65') 

(L w /a) 4a + (L w /a) 2a ’ 1 j 


6.2 Conductances under RG flows 

We will now use our understanding of the RG flows to 
study the conductances of a three-wire system as func¬ 
tions of physical parameters such as the wire lengths and 
the temperature, for the case of repulsive interactions. 
In particular, we will study how the conductances scale 
with various lengths when we approach the completely 
stable fixed point. 

Using the procedure described in Eqs. HELP and (ED), 
we first find the values of /3 for different flow directions 
given by the perturbation Si around a fixed point So- 
The stable fixed point that we are interested in has only 
the r ee jj ’s and rhh,jj ’s (namely, the diagonal elements of 
So) being unimodular numbers and all the other elements 
being zero. For a flow direction Si in which only the 
phases of the diagonal elements are changed, we find that 
de/dl = 0. This implies that /? = 0 and the RG flow is 
marginal. This is expected since the RG fixed points 
are invariant under multiplication by a diagonal unitary 
matrix as discussed in Eq. (15(71) ; hence there is no RG flow 
in those directions. Next, we look at the RG flows when 
only the r e h,jj’s and r^ejj’s are perturbed from zero. 
For this perturbation, we find that de/dl = — 2ae; hence 
/3 = —2a, and the RG flow in this direction is irrelevant. 
We choose Si where either (i) only the t e e,y’s and thh,ij 's 
are non-zero, or (ii) only the t e h.ij's and th e ,ji s are non¬ 
zero. For both these cases we find that de/dl = —ae, 
which means that /3 = —a and the RG flows in these 
directions is also irrelevant. The RG flows near the fixed 
point So are therefore either marginal or irrelevant in all 
directions. 

Now we will discuss how the Cooper conductance Gc 
and the normal conductances Gn 2 and Gn 3 scale under 
the RG flows (we are assuming that an electron is in¬ 
cident from the NM lead 1). Physically there are three 
length scales in the problem and we have to stop the 
RG flows when we reach the smallest of the three scales. 
One length scale is 77 = Hvf/A which is associated with 
the SC gap, another is the wire length L w (we will take 
the lengths of all the three wires to be of the order of 
L w ), and the third scale is the thermal coherence length 
Lt = Hvf/ (k B T) if the system is at a temperature T 85 ' 86 . 
(We assume that all these length scales are much larger 
than the short distance scale a). We will now consider 
different regimes of these length scales. 

(i) We first consider the case in which L w is finite and 
smaller than 77 , while T —> 0 so that Lt —> oo. Then 


where Ci, c 2 are some constants which depend on how 
far from the stable fixed point we are at the length scale 
a where the RG flows begin. If L w a, the second term 
in Eq. (ESI) dominates over the first term. We therefore 
obtain 


Gc 


1 

{L w /o) 2a ' 


( 66 ) 


Similarly we find that Gn 2 , Gn3 ~ \t e e \ 2 — \the \ 2 scale 
as 


Gn2, Gn3 


1 

{L w /a ) 2a " 


(67) 


(ii) Next we consider the case in which T is finite such 
that Lt = hvp / (kpT) is smaller than 77 , and L w —> 00 . 
Then the RG flows stop at the length scale Lt- Hence 
the elements of S must be evaluated at l = ln(Lp/a) = 
ln(vF/(k B Ta)). Hence 

G c ~ diksTa/Hyp) 401 + c 2 {k B Ta/hv F ) 2a - (68) 


In the limit ksTa/h-vp 1, the second term dominates 
over the first and we get 

G c ~ ( k B Ta/hv F ) 2a - (69) 

Similarly, we find 

Gn 2 , Gn3 ~ ( k B Ta/hvF) 2a ■ (70) 


(iii) Finally we consider the case where 77 = Hvf /A is 
smaller than both L w and Lp. Then the RG flows must 
be stopped at the length scale 77 since our RG equations 
were derived under the condition that the energy scale is 
much larger then A. We then obtain expressions for the 
conductances which are similar to Eqs. (lEH) and (1671) but 
with Lw replaced by 77 . 

To conclude, we see that the conductances Gc and 
Gnj approach zero as powers of the smallest length of 
the system which may be L Wl Lp or 77 . The exponents 
of the power laws can give an estimate of the strength of 
the interactions in the wires. 


7 Experimental realization of systems 
with different signs of A j 

We will now discuss how it may be possible to experi¬ 
mentally realize a system of three SC wires with the same 
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or different signs of the p-wave pairings A j. We consider 
the system studied in Ref. El This consists of a wire 
with a Rashba spin-orbit coupling of the form ±app r a x , 
where p r is the momentum along the wire and a x is a 
Pauli spin matrix. This form can arise as follows. Let us 
take the coordinate in the wire to increase along an arbi¬ 
trary direction f lying in the x — y plane. If the Rashba 
term is apfi ■ a x p, and h points in the z direction, then 
the Rashba term will be app r a x if f = y and —anp r a x 
if f = —y. 

Next, we place the wire in a magnetic field in the z 
direction which is perpendicular to the Rashba term (and 
has a Zeeman coupling Az to the spin of the electrons) 
and in proximity to a bulk s-wave SC with pairing Ag- 
The complete Hamiltonian isi^ 


H 



y)8 a p ± a R p r cjap 


A z o z a p 




+ l - J dr +H.c.}, 


(71) 


where \I , Q , is the annihilation operator for an electron with 
spin a, and the ± sign of the Rashba term depends on 
whether f = ±y. Ref. |H then shows that for a certain 
range of the parameters, this system is equivalent to a 
spinless p-wave SC of the form that we have studied in 
this paper, with the p-wave pairing term being given by 
~i{A/kF){c i d x d J r d)d x c) (see Eq. (JT]) ), where 


A 


ankpAg 

Az 


(72) 


Now consider a case in which the coordinates are given 
by f = y in all the three wires; see Fig. [5] (a). (The region 
where the three SC wires meet is an extended vertical 
region on the left side of the figure). Then the Rashba 
term and hence Aj will have the same sign in all the SC 
wires. On the other hand, suppose that one of the SC 
wires runs in a direction which is opposite to the other 
two wires as shown in Fig. [5] (b). (The three SC wires 
now meet along an extended vertical region in the middle 
of the figure). Now f = — y in that wire while f = y in 
the other two wires. Then Eq. C72D shows that Aj will 
have one sign in that wire and the opposite sign in the 
other two wires. 

[We would like to note that if the angle between the 
wires is different from zero or 7r, the situation will be 
more complicated because the Rashba term apn ■ a x p 
will no longer be proportional to the same cr matrix in 
the different wires. Hence the effective p-wave pairings 
in the different wires will not be related simply by sign 
changes in Aj}. 

It is clear that the vertical region where the three SC 
wires meet is likely to cause some scattering of the elec¬ 
trons; we have modeled this scattering in the earlier sec¬ 
tions using the matrix M. Finally, the ends of the SC 
are connected to NM leads through tunnel barriers. As 
discussed in Sect. 2, these barriers can be characterized 
by their strength A. 


In Sects. 2 - 4, we discussed a conductance Gc in which 
pairs of electrons can appear in (or disappear from) the 
SC regions. At a microscopic level we can understand 
these processes as occurring due to a Cooper pair going 
from the bulk s-wave SC to one of the p-wave SC wires 
(or vice versa). Finally we assume that the s-wave SC 
is grounded, and the three NM leads and the s-wave SC 
form a closed electrical circuit so that we can measure 
the conductances Gc and Gjvj- 

A 


A 


A 


(a) 


(b) 

FIG. 5: Three SC wires with (a) the sign of A j being the 
same in all the wires, (b) the sign of Aj being different in 
the wire on the left compared to the other two wires. The 
SC wires are shown in a lighter shade while the NM leads are 
shown in a darker shade. 


8 Conclusions 

In this paper, we have studied the conductances of a 
system consisting of three p-wave superconductors, each 
connected to a normal metal lead (these are labeled as 
SCI, SC2, SC3 and NM1, NM2, NM3 respectively) in two 
different regimes of the energy. In Sects. 2-4, we have 
studied the Majorana and related sub-gap modes and 
their effects on the sub-gap conductances. Here we have 
considered two cases: (i) when the p-wave pair potentials 
Aj have the same sign in all the SCs, and (ii) when one of 
the Aj has a different sign from the other two. In Sects. 
5 - 6, we have studied the effect of interactions between 
the electrons on the conductances at energies lying far 
from the SC gap. 

To study the sub-gap modes, we have used a continuum 
model and current conservation to derive the boundary 
conditions at the junctions between the NMs and the 
SCs. Then we have found the boundary condition at the 
junction of the three SCs; this condition is encoded by 
a Hermitian matrix M. Using these conditions, we have 
numerically studied two conductances, the Cooper pair 
conductance Gc (from NM1 to the SCs) and the normal 
conductance Gn 2 (from NM1 to NM2) when we send in 
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an electron from NM1 with an energy E; we have taken E 
to he inside the superconducting gap so as to specifically 
probe the sub-gap inodes which reduce to zero energy 
Majorana modes in the limit of long wire lengths. 

We have looked at three different regimes of the lengths 
of the SCs with respect to rj which is the length scale 
associated with the SC gap. We find a rich pattern of 
the conductances as functions of the SC lengths and E. 

We first consider the case when all the Ay’s have the 
same sign. In the length regime L ~ r) we find six sub¬ 
gap modes; these are not at zero energy because they 
hybridize with each other due to the finite lengths of the 
SCs. Among the six modes, three sit near the SC-NM 
junctions, and the other three sit at the junction of the 
three SCs. All these modes have a significant effect on 
the conductances and the latter shows peaks exactly at 
the energies of the sub-gap modes. 

For the length regime L ~ 2rj the hybridization among 
the sub-gap modes is less than in the previous case; hence 
the energy splitting of the sub-gap modes is also less. 
Even if we cannot clearly see six separate sub-gap modes, 
we can still observe several modes at non-zero energies. 

In the third regime, when L>j), we find that the sub¬ 
gap modes are almost decoupled from each other and 
hence lie at zero energy. The Cooper conductance Gc 
almost approaches its highest value of 2 e 2 /h. 

We compare our conductance results with the sub-gap 
energies that we get for a box made of three SC wires 
with hard walls (namely, without the NM leads). We 
generate plots for the energies of the sub-gap inodes in 
all the three length regimes as discussed earlier. In the 
first length regime, instead of six sub-gap modes we see 
only four clearly; we believe this is due to the fact that 
two of the modes are very close to E/A = ±1 and hence 
lie beyond our resolution. In the other two regimes, the 
plots of the energies of sub-gap modes almost exactly 
match with the conductance plots for the full system as 
we can see in Fig. [2] 

In obtaining all these results, we have made a partic¬ 
ular choice of the matrix M which defines the boundary 
condition at the junction of three SCs, namely, we have 
taken all the diagonal and off-diagonal elements of M to 
be equal to 1. If we take the elements to be different, the 
results may be somewhat different from what we get. 

With the same choice of M, we have calculated the 
conductances taking the A in SCI to have a different 
sign compared to the A’s in SC2 and SC3. The results 
are significantly different from the case of all Ay’s having 
the same sign. In the length regime L ~ rj we find four 
sub-gap modes at different energies, instead of six. Now 
there is only one sub-gap mode sitting at the junction 
of three SCs. As we increase L/rj, the energy splitting 
between the sub-gap modes decreases; when L > ?; we 
find that all the modes he at zero energy. Once again we 
have compared our conductance results with the energies 
of the sub-gap modes in a SC box with the A in SCI 
having an opposite sign to SC2 and SC3. The results 
match very well as we can see in Fig. [3] 


We have then presented some analytical and symmetry 
arguments to explain the presence of multiple zero energy 
Majorana modes at a junction of three long SCs. Using 
the boundary condition at the junction we find that the 
number of independent variables in the problem directly 
corresponds to the number of Majorana modes at the 
junction. We have also used an “effective time reversal 
symmetry” to argue that the number of Majorana modes 
at a junction of three Kitaev chains can be either one or 
three. 

Next, we have studied the effect of interactions between 
the electrons in order to understand the conductances 
at energies far from the SC gap. (We emphasize that 
the energy range we have considered is different from 
earlier work where RG equations have been studied at 
energies within or close to the SC gap, where the various 
reflection and transmission amplitudes vary rapidly with 
the energ y 16 ! 90 ). We have derived the RG equations for 
a general scattering matrix S which governs a system of 
several NM wires which meet at a junction with a SC 
region. The RG flows of the elements of S are entirely 
a result of the interactions in the NM wires; there is no 
flow if all the interaction parameters ay are equal to zero. 
We have found stable and unstable fixed points for two- 
wire and three-wire junctions when the interactions are 
repulsive, i.e., ay > 0. A completely stable fixed point is 
one where the RG flows take us back to the fixed point for 
a small deviation from it in any direction. We find that 
at the completely stable fixed point, the magnitudes of all 
the normal reflection probabilities are equal to 1, while 
at the completely unstable fixed points all the Andreev 
reflection probabilities are equal to 1. We have discussed 
the consequences of the RG flows on the conductances of 
the system, assuming all the interaction parameters to be 
equal to a. A stability analysis near the stable fixed point 
shows that when the wire lengths L w are finite and the 
temperature T — > 0, the Cooper conductance Gc and 
the normal conductances Gn 2 , Gn 3 all scale as 1 / L 2 °-. 
In the other limit where T is finite (but small compared to 
the band width of the system), and L w —> oo, Gc, GN 2 
and Gn 3 scale as T 2a . In general, for an electron incident 
from NM lead j, the conductance measured on NM lead 
j will scale as 1 /Li/// or T 2a depending on which length 
scale is smaller. Thus a measurement of the conductances 
can provide valuable information about the strength of 
the interactions between the electrons. 

We can summarize our most important results as fol¬ 
lows. 

(i) We have shown that a system of three p -wave SC 
wires with NM leads has multiple sub-gap modes and all 
these modes contribute to peaks in the conductance; the 
positions of the peaks exactly match the energies of the 
sub-gap modes. We thus have a novel system with a large 
number of sub-gap conductance peaks which vary with 
the applied bias and wire lengths in an interesting way. 

(ii) To study the conductances far from the SC gap, 
we have used an RG method which directly uses the 
fermionic language to look at the effect of interactions 
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on the scattering matrix. (This is in contrast to ear¬ 
lier RG studies of p -wave superconductors which used 
bosonization and studied the RG flows of parameters in 
the bosonic Hamiltonian rather than the scattering ma¬ 
trix) . We have shown that interactions make the scatter¬ 
ing matrix flow to certain stable fixed points at large wire 
lengths and low temperatures. The fixed points are sym¬ 
metric under permutations of the wires and under time 
reversal even if the system does not have these symme¬ 
tries at the microscopic length scale. Near the stable 
fixed point, we have shown that the conductances scale 
as a power of the wire length or the temperature, and 
the power can give us an estimate of the strength of the 
interactions. 

Putting the above results together, we get a complete 
picture of the conductances of a system of three p -wave 
SC wires both inside and far outside the SC gap. 

Finally, we have discussed how our model can be ex¬ 
perimentally implemented by taking three wires with 
Rashba spin-orbit coupling, applying a Zeeman field per¬ 
pendicular to the direction of the Rashba field, and plac¬ 
ing the system in proximity to an s-wave superconductor. 
In a particular range of parameters, each wire effectively 


becomes a p -wave SC. The cases where the three SC wires 
have the same sign of the p -wave pairing or one of them 
has the opposite sign of the pairing can both be imple¬ 
mented, simply by rearranging the orientation of one of 
the wires. As we have shown in the earlier sections, the 
number of sub-gap modes is different in the two cases, 
and this leads to a significant difference in the pattern of 
conductance peaks as a function of the energy and wire 
lengths. In future, we can generalize our studies of a 
single three-wire junction to more complicated systems 
such as networks of wires forming a lattic o 11 ! 66 , where 
each vertex of the lattice can host one or more sub-gap 
modes. 
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